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I. INTRODUCTION 


After the discovery of the Higgs (BEH) boson JT| 2], all fundamental particles of the 
Standard Model (SM) are finally found, and now even passionate adepts of the SM should 
look for physics beyond it. The pattern of particles we have is rather asymmetric: there are 
twelve vector bosons, many leptons and quarks with spin 1/2 and only one scalar particle 
h with mass 125 GeV. Of course, there is only one particle with spin 2 as well, a graviton. 
However, unlike the spin 2 case, there are no fundamental principle according to which 
there should exist only one fundamental scalar particle. That is why it is quite probable 
that there are other still undiscovered fundamental scalar particles in Nature. The purpose 
of the present paper is to consider the simplest extension of the SM by adding one real scalar 
field to it. Such an extension of the SM attracts considerable attention: relevant references 
can be found in recent papers M. Extra singlet can provide first order electroweak phase 
transition needed for electroweak baryogenesis. It can act as a particle which connects SM 
particles to Dark Matter. Not going into these (very interesting) applications, we will study 
the degree of enhancement of double higgs production at LHC due to an extra singlet. To do 
this we should analyze bounds on the mass of the additional scalar particle and its mixing 
with isodoublet state. 


An enhancement of hh production occurs due to the mixing of the SM isodoublet with 
additional scalar field which is proportional to the vacuum expectation value (vev) of this 
field. Thus isosinglet is singled out: its vev does not violate custodial symmetry and can be 
large. For higher representations special care is needed; see paper [7] where an introduction 
of isotriplet (s) in the SM is discussed. 


The paper is organized as follows: in Section [TT] we describe the model and find the 


physical states. In Section [ITT] we get bounds on the model parameters of the scalar sector 
from the experimental data on h production and decays and from precision measurements 
of Z- and IL-boson parameters and /-quark and h masses. In Section [TV] we discuss double 
h production at LHC Run 2. In Appendix [A] quantitative description of single and double 
higgs production at LHC is presented. 
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II. THE MODEL 


Adding to the SM a real field A", we take the scalar fields potential in the following form: 

V($,X) = + ^AA 2 + (1) 

where $ is an isodoublet[j] Terms proportional to A 3 , X 4 and <f)f$A 2 are omitted despite 
that they are allowed by the demand of renormalizability: we always may assume that they 
are multiplied by small coupling constants. Two combinations of the parameters entering (|TJ) 
are known experimentally: it is the mass of one of the two scalar states, h, which equals 
125 GeV and the isodoublet expectation value = 246 GeV. The two remaining combina¬ 
tions are determined by the mass of the second scalar, H (we take mn > though this is 
not obligatory), and the angle a which describes singlet-doublet admixture: 


Substituting in (|TJ) 


h = (f) cos a + x sin a, 
H = —(f) sin a + x cos a , 

$ = 


,72 (vv + t + iv), 

at the minimum of the potential we get: 


(f> — h cos a — H sin a, 
X = h sin a + H cos a. 

X = v x + X, 


{ An| + 2 jjLv x = ml, 
2 m\v x + fJ-v | = 0, 


so /r is negative. For the mass matrix using (|dj) we get: 

Vchd, !4v \ / Av| 


M = 


'(jxj) V(j>x 

V<t>x V XXi 


liv® m x 


where V^ x = ■ ■ ■ Eigenvalues of ([5]) determine masses of scalar particles: 

™1,h = + \n 2 x T 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


1 We are grateful to J. M. Frere who brought to our attention that similar model was considered long ago 
in 0. 
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FIG. 1: Dependencies of the model parameters on the mixing angle for = 300 GeV. 


where ” corresponds to m,h and “+”—to rn h■ Eigenfunctions are determined by the 
mixing angle a\ 


sin 2a 


-2 nv® m 2 h - Aw| 

—5 -tan a— - 

m z H — ml 


(7) 


Equations ([7]) determine /i and A for the given mixing angle a, while equations © 
determine mx for given a as well. Finally, equations (Ul) determine the values of and Vx- 
Fig.0 demonstrates the dependencies just described for m H = 300 GeV. 
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III. BOUNDS FROM h PRODUCTION AT LHC AND ELECTROWEAK PRECI¬ 
SION OBSERVABLES 


ATLAS and CMS collaborations had detected h production and decays in the reactions 


PP -» h -» fi, 


( 8 ) 


where fi,i = 1, 2,..., 5 designate the so-called “Big five” final state channels: WW*, ZZ*, 
77 , tt, bb. Cross sections of reactions ([ 8 ]) are equal to the higgs production cross section 
times branching ratio of the corresponding decay channel. Quantities p t are introduced 
according to the following definition: 


— a pp^h ' r^y./Tfe 

(<w r h ^ fi /r h ) su 1 J 

According to ATLAS and CMS results, all /q are compatible with one within experimental 
and theoretical accuracy. It means that no New Physics are up to now observed in h 
production and decays. 

In the model with an extra isosinglet, production and decay probabilities of h equal that 
in the SM multiplied by a factor cos 2 a, that is why we have: 


Pi = cos 2 a , 


( 10 ) 


and existing bounds on p t are translated into bounds on the mixing angle a. Taking into 
account all measured production and decay channels, for the average values experimentalists 
obtain PEI!: 

ATLAS: p = 1.30±g;i|, (11) 

CMS: p = 1.001q!i3 [±0.09(stat.)^oor(theor.) ± 0.07(syst.)] (12) 

Let us stress that the theoretical uncertainty in the calculation of pp —» h production cross 
section at LHC does not allow to reduce substantially the uncertainty in the value of p. 
Bounds from electroweak precision observables (EWPO) are not affected by this particular 
uncertainty. 

We fit experimental data with the help of LEPTOP program m using mu = 125.14 GeV. 
The result of the SM fit which accounts the h mass measurement is shown in Table [Tj Quality 
of the fit is characterised by the x 2 value 

X 2 Au.o.f. = 19.6/13. (13) 
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TABLE I: EWPO fit of the Standard Model 


Observable 

Experimental value 

Standard Model 

Pull 

r z , GeV 

2.4952(23) 

2.4966(14) 

-0.5895 

a h , nb 

41.541(37) 

41.475(14) 

1.7746 

Ri 

20.771(25) 

20.744(18) 

1.0831 

A 1 

^FB 

0.0171(10) 

0.0165(2) 

0.6572 

A r 

0.1439(43) 

0.1484(7) 

-1.0452 

Rb 

0.2163(7) 

0.2158(0) 

0.7699 

R c 

0.1721(30) 

0.1722(0) 

-0.0277 

A b 

^FB 

0.0992(16) 

0.1040(5) 

-3.0303 

A c 

^FB 

0.0707(35) 

0.0744(4) 

-1.0565 

s 2 (Qfb) 

0.2324(12) 

0.2313(1) 

0.8771 

Alr 

0.1514(22) 

0.1484(7) 

1.3822 

Ab 

0.923(20) 

0.9349(1) 

-0.5941 

A c 

0.670(27) 

0.6685(3) 

0.0567 

M w , GeV 

80.3846(146) 

80.3725(67) 

0.8322 

rrit, GeV 

173.24(95) 

174.32(89) 

-1.1370 

1/a 

128.954(48) 

129.023(37) 

-1.4378 


Higgs boson contributions to clectroweak observables at one loop are described in LEP- 
TOP by functions H t (h) = H i {m L h jrn 2 7 ). In the case of an extra singlet the following substi¬ 
tution should be performed: 

Hi(h) —> cos 2 a Hi(h ) + sin 2 a H^H), H = m 2 H /m 2 z . (14) 

The same substitution should be made for the functions S^Viitjh), t = which 

describe two loops radiative corrections enhanced as m\. In two loops quadratic dependence 
on higgs mass appears which is described by functions <f 5 V). Calculations of these corrections 
in the case of an extra singlet higgs is not easy. An approximate upper bound has been 
estimated by assuming that 

5 5 Vi(H) < <^((1000 GeV) 2 /m|) « 100 S 5 Vi(h) for m H < 1000 GeV. (15) 
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Comparison of two calculations, one with 8 ^Vi{h) = cos 2 a 8 ^Vi{h), and the other with 


85 Vi(h) = cos 2 a 8 ^Vi{h) + 100 • sin 2 a S^V^h), (16) 

showed that the correction to the values of sin a in Fig. [ 2 ] is less than 10 -3 . 

Bounds from EWPO on the singlet model parameters are presented in Fig. 2a % 2 mini¬ 
mum is reached at sin a = 0, = 150 GeV, which is the minimum value allowed for m# 

in the fit. Experimental data are avoiding heavy higgs. The value of y 2 at the minimum 


coincides with the SM result (13). Lines of constant y 2 correpospond to Ay 2 = 1,4, 9,.... 
Probabilities that (sin a, run) values are below these lines are 39%, 86%, 98.9%, ... |^] 
Bounds accounting for both EWPO and direct h production data ([TTJ) , ( [l2| are shown 
in Fig. 2b We see that for heavy H bounds from EWPO dominate, while for light H 
measurement of fj, is more important. 


IV. h, H AND hh PRODUCTION AT LHC 


The main purpose of this section is to find what enhancement of double higgs production 
cross section is possible with enlarged higgs sector. Let us remind that in the SM double h 
production cross section is very small. According to the recent result [TUI, at a/s = 14 TeV 
a NNLO (pp —y hh) = 40 fb with a 10 A 15% accuracy. We will demonstrate that enlarged 
higgs sector allows to strongly enhance double h production. 

The cross section of H production at LHC equals that for the SM higgs production (for 
(mh )sm = m#) multiplied by sin 2 a. Cross section of the SM higgs production at NNLO we 
take from Table 3 of [13]. In order to obtain cross section of resonant hh production in H 
decays we should multiply cross section of H production by Bt(H —>■ hh). 

Let us consider H decays. Decays to hh, W + W ~, ZZ and tt dominate. For the Hhh 
coupling we obtain: 


A Cnhh — 


cos 2 a sin a — ^ cosa/l — 3 sin 2 a.) 
2 ml + m 2 H 2 


Hh 2 


2u$ 


sin a cos 2 a Hh" 


(17) 


= gHhhHh 2 , 


2 Let us note that if a subset of experimental data from Table ||] is fitted, then allowed domains of the 
(sin a, mu) values will be larger than those presented in Fig. |2a| Here we disagree with the statement 
made in [4] that the fit of only one observable ( mw ) allows to set the strongest constraint on (sin a. niff ). 
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(a) Bounds from electroweak precision observables. 
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(b) Bounds from both electroweak precision observables and signal strength 
measurements The dashed line corresponds to Ax 2 = 5.99; the probability that 

numerical values of {run, sin a) are below it equals 95% (compare with Ref. [6], eq. (23)). 


FIG. 2: Bounds on the singlet model parameters. 










sm a 

(a) Decay widths. 


sm a 

(b) Branching ratios. 


FIG. 3: Decay widths and branching ratios of the heavy higgs boson for = 300 GeV. 
thus 
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9Hhh 




2 mv 


8nrriH V V 771 h 
Decays to W + W~, ZZ, tt occur through isodoublet admixture in H : 


(18) 


a c = 2 Hk 


sm a 


HW + W- + ^ — - ur72 ' mt 


= gnwwHW+W- 


sin a HZ 1 H- 1 sin a Htt 

w<f, 

OhzzHZ 2 + g m iHti } 


(19) 


thus 


r h^w+w — 


H^ZZ — 


9hww 7JI h 

9hzz 7TI h 

12871771 ^ L 


1-4- 


, m 


w 


+ 12 - 


m 


w 


m 


H 


m 


H 


fW) 


1-4- 


m 


+ 12 - 


rrir 


m 


H 


m 


H 


2 m z \ 2 

m H J 


r 


H^tt 


3g 2 m?riH_ 

8 tt 


1 - 


2 m t 
m H 


( 20 ) 

( 21 ) 

( 22 ) 


The dependence of the widths and branching ratios of H decays on mixing angle a for 
rriH — 300 GeV are shown in Figure |3j 

For the cross section of the reaction pp —> H —> hh we have: 


a(pp —> H —> hh) = o(pp —> h ) SM • sin 2 a ■ Br (H —>■ hh), 


(23) 
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FIG. 4: Contour plot of a(pp —> H —> hh) for y/s = 14 TeV. 
In this figure we neglect small effects of H —> hh*. 
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FIG. 5: Contour plot of R = . 

n (a(pp^h)Br(h^ZZ)) SM 

In the calculation of R we assume m# > 2ro ft . 
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the lines of constant cross section are shown in Fig. [4] (compare to Fig. 4 from [6j). H —» ZZ 
decay can be used in order to find 77; its cross section divided by that for the SM higgs 
boson with (rrih) sm = m H is 


a(pp ->■ 77) • Br(77 -> ZZ) 
(cr(pp h ) ■ Br (h -)■ ZZ)) SM 


sin 4 a 

• 2 , r(H->hh) ’ 

sm a H— L 

1 SM 


(24) 


Contour plot of R is presented in Fig. [HJ Let us note that R does not depend on 


s. 


V. CONCLUSIONS 

In the models with extended higgs sector strong resonant enhancement of double higgs 
production is possible which makes the search of pp —y hh reaction at Run 2 LHC especially 
interesting. According to Fig. [4] cross section of pp —y H —> hh reaction can be as large as 
0.5 pb, ten times larger than the SM value. 

The search for H boson can also go in the same way as it was for the heavy SM boson 
h. Probability of H observation diminishes compared to that of h because of a) suppression 
of H production cross section by the factor sin 2 ex < 0.2; b) suppression of Br(iJ —» ZZ) 
because of additional H —* hh decay mode. Taking these two factors into account, we get 
about factor 10 suppression of pp —>■ H —>■ ZZ process probability compared to that for the 
SM higgs boson (see Fig. [5]). 

Results for the search of higgs-like boson in ZZ decay mode can be found in [20], Figure 5. 
Comparing it with our Fig. [5] we observe that experimental data start to be sensitive to the 
singlet model expectation for maximally allowed values of the mixing angle a. 

After the first version of this paper was published in arXiv, we got a number of emails 
providing us with references to related research [2TJ. 

S. G., M. V. and E. Zh. are partially supported under the grants RFBR No. 14-02-00995 
and NSh-3830.2014.2. S G. and E. Zh. are also supported by MK-4234.2015.2. In addition, 
S. G is supported by Dynasty Foundation and by the Russian Federation Government under 
grant No. 11.G34.31.0047. 
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Appendix A: Higgs production in effective Lagrangian approach 


Simple analythical formulas which qualititavely describe single and double higgs produc¬ 
tion in the SM are presented in this section. Let us start with single higgs production in 
gluon fusion. In the limit rrp,, <C 2m*, the amplitude of gg —y h transition is determined by 
the top quark contribution into the QCD Gell-Mann-Low function: 


A£ = 


(y. s 

1271 


In ( 1 + — ) G 


M = -^G'GLh, 


fils') 


6 nV§ 


Ills'-* fils' 


leading to the well-known result for the production cross section: 


0% T 0 X , \ 
a "^ h 5767ru 2<5(r To) ’ 


(Al) 


(A2) 


Here r = s/s and To = rri^/s; s = (j)\ + P 2) 2 is the invariant mass of colliding protons, 
s = X\X 2 s = ts is the invariant mass of colliding gluons. Integrating over gluons distribution 
in a proton, we get: 

1 1 

a pp ^h = J dx x I dx 2 g{x l )g(x 2 )(Tgg^ h . (A3) 

r 0 T 0 /x 1 

Changing the variables from aq, x 2 to r, y according to the following definitions: aq = a /7e y , 


x 2 = s/re y , and substituting (A2) into (A3), we obtain: 

-ln^/ro 

aX 1 f 
p p ^ h 5767 tu | s J 

In Vtd 


g(^)g(^- y )dy = dL 


5767 tu| ds ’ 


(A4) 


where the so-called gluon-gluon luminosity is given by the integral over gluon distributions: 

-ln^ro 

= - [ g(V^o eV )9{V^e~ y )dy. (A5) 

nl S J 

h In 


dL 

ds 


A number of PDFs parametrizations exist in the literature; their results for (A5) at 
a/s = 7, 8, 14 and 100 TeV and m\ = (125 GeV) 2 coincide within several percents. Finite 
value of m t = 172 GeV should be taken into account by multiplication of the leading order 


result for the amplitude M (Al) by a factor 


F = ^[(1 - /V + 1], 


where /3 = 


2 m/ 


, and x = arctan 


for (3 > 1, x = - ( 7r + % In 


(A6) 

1 + a/1 — fd 


mi, J \J3 ~ 1 2 

for /3 < 1 [Qj (note that lim F — 1). This adjustment leads to 6% enlargement of <x gg ^,h 

m*—>-oo 
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FIG. 6: Leading-order diagrams for the double higgs production at LHC. 


compared to m t —> oo value; however taking into account b and c quark contributions results 
in 6% overall reduction. 

Applying all these factors and using PDFs from S2. we obtain numbers presented in 
Table [IlJ To calculate cr NNLO from a LO we use K -factor from [13]: K m 2.5 for y/s = 7 and 
8 TeV, and K « 2 for yfs = 14 TeV. For y/s = 100 TeV K « 1.5 (A. Djouadi, private 
communication). Let us stress that according to [T3] , accuracy of the calculated value of 
AA? is at the level of ±(10 ± 17)% which makes hopes of reducing uncertainty in /i ?; (and 
p) below 10% elusive. In the case of an extra singlet, h and H production cross sections 
equal the SM one multiplied by cos 2 a and sin 2 a respectively. 

Let us turn now to double h production at pp collision in the SM. At the leading order it 
is described by the two diagrams shown in Fig. [6} According to equations (4) and (11) and 
Table 1 from im the cross section of the double production of the 125 GeV h at 14 TeV 
LHC in the leading order equals: 

a LO (pp -y hh) = 144.6 • (0.169 2 + 0.457 2 - 1.79 ■ 0.457 • 0.169) fb = 14 fb, (A7) 


where the first term in parentheses originates from the square of the triangle diagram, the 
second—from the square of the box diagram, while the last one is their interference, which 
diminishes the cross section. 


In order to understand result (A7) let us proceed in the following way. In the limit 
s 4m 2 the triangle gg —)• h and box gg —y hh amplitudes can be directly extracted from 


lagrangian (Al), expanding it over h /tq,: 

AC = In 


12tt 


h 

WrJ, 


a* 


1 + -'±-= 12 , 


A _ 1^! l r 2 

WrJ, 2 V 2 > ^ 


(AS) 


where the first term corresponds to the diagram shown in Fig. |6a} while the second term 
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TABLE II: Data relevant for the SM higgs boson production at LHC. The difference 
between the numbers in Tables He and |IId is due to poor accuracy of K -factors presented 


m 


(a) §, IQ" 3 GeV 2 . 


\ VS 

m H \ 

7 TeV 

8 TeV 

14 TeV 

100 TeV 

125 GeV 

6.41 

8.30 

22.9 

451 

300 GeV 

0.147 

0.205 

0.737 

25.1 


(b) a LO (pp —> h ), pb. 


\ VS 

m H \ 

7 TeV 

8 TeV 

14 TeV 

100 TeV 

125 GeV 

5.52 

7.16 

19.8 

389 

300 GeV 

0.936 

1.31 

4.69 

160 


(c) cr NNLO (pp —> h ), pb. 


\vs 

rn H 

7 TeV 

8 TeV 

14 TeV 

100 TeV 

125 GeV 

13.8 

17.9 

39.6 

583 

300 GeV 

2.34 

3.27 

9.37 

239 


(d) cr NNLO (pp —> h), pb, from Tables 1, 3 of [T5] . 


\ VS 

m H \ 

7 TeV 

8 TeV 

14 TeV 

100 TeV 

125 GeV 

15.31 

N/A 

49.97 

N/A 

300 GeV 

2.42 

N/A 

11.07 

N/A 


describes the diagram shown in Fig. [6I3J Triple higgs coupling is given by the following term 
in the SM lagrangian: 

AC = (A9) 

which leads to A hhh = 3 m\/v§. Hence, for the sum of the triangle and the box diagrams at 


14 



































s <C Am 2 we get 


M = 


CXa 


Q'KVs 


m 


s — mi 


3 n 


1 s^i 2 


(A10) 


which equals zero at threshold when s = (2ro/,) 2 [18, IS]. For the cross section we get 


Vgg^fhh \ 


a 2 G 2 F s 


K<4m * 576(2 tt) 3 


1 - 


3 m 


2 F 

h 


s — m 


i - 


(2 m h y 


(All) 


(see Eq. 13 from [IS]). 

In the high-energy limit s Am 2 box diagram dominates and the cross section behaves 


as: 


Ogg^thh \ 

s>4m? 


= A* 


167t 3 s V 


— I \H 


(2 m h y 


(A12) 


Normalization constant A is determined by the condition that at s — Am 2 expressions (All) 


and (A12) are equal: 


A = - 
6 


3 m 2 h 


Am? — m 


(A13) 


Finally, for the cross section of double h production in the SM we obtain the following 
approximate expression: 


&pp—>hh 


dS (Tgg^fhk (^) 1 


(2 m h ) 2 

— In Jr 


dL _ 1 
ds s 


g(\ffe v )g(\/Te y )dy, 


(A14) 


(A15) 


In y/r 


where Equations (All)—(A13) should be substituted in (A14) and r = s/s, s being the hh 
invariant mass. The differential cross section is shown in Fig. [7], while for the total cross 
section for hh production in the SM we get a{pp —>• hh) = 4 fb at a/s = 14 TeV, 3.5 times 


smaller than the explicit leading order result (A7). 
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Appendix B: Colored figures 
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Bounds from electroweak precision observables. 
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2 b 


Bounds from both electroweak precision observables and signal strength measure¬ 
ments (11), (12). The dashed line corresponds to Ay 2 = 5.99; the probability that numerical 
values of (m#, sin a) are below it equals 95% (compare with Ref. [6], eq. (23)). 
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Bounds on the singlet model parameters. 
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rn h , GeV 

FIG. [lj Contour plot of a(pp —> H —» hh) for a/s = 14 TeV. In this figure we neglect small 
effects of H —» hh*. 




FIG. [5I Contour plot of i? = ^ . In the calculation of f? we assume 

rn h > 2 m/,,. 
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